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Abstract. We study the impact of an intracluster magnetic field on the main structural properties of clusters 
and groups of galaxies: the radial density and entropy profiles, the entropy - temperature relation and the X-ray 
luminosity - temperature relation for groups and clusters of galaxies. To this aim, we develop a description of 
' the intra-cluster gas based on the Hydrostatic Equilibrium condition and on the Magnetic Virial Theorem in the 

I/") , presence of a radial distribution of the magnetic field B(r) — B,(p g (r)) a , with a w 0.9, as the one indicated by 

observations and numerical simulation. Our analysis shows that such a description is able to provide, at once, a 
possible explanation of three problematic aspects of the cluster structure: i) the flattening of the entropy profile 
in the cluster center ii) the flatness of the S — T relation iii) the increasing steepening of the Lx — T relation 
from the cluster scale towards the group scales. The available entropy and X-ray luminosity data indicate that 
an increase of the magnetic field B, ~ y0.5±o.i j g re q U j re( j ^o reproduce at the same time both the S — T 
and the Lx — T relations. It follows that a consistent description of the magnetized ICM can provide a simple 
CIh' explanation of several (or of all) of these still open problems, and thus weakens the need for the inclusion of 

other non-gravitational effects which have been proposed so far for the explanation of some of these features. This 
(initial, but not conclusive) analysis can be regarded as a starting point for a more refined analytical exploration 
of the physics of the magnetized intra-cluster medium, and it provides testable predictions that can be proven or 
disproven with the next coming sensitive observations of groups and clusters in the X-ray band and in the radio 
frequency band. 
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1. Introduction perature, entropy, X-ray luminosity, mass) to the tem- 
perature of the IC gas (i.e., a measure of their poten- 

It has been well known since the early 1990s that the intra- ^ ^ &nd q{ their equilibrium stage) show that thc 

cluster (IC) gas properties deviate from the self-similar glob&1 JC gag properties nke the X -ray luminosity L x 
predictions of the standard CDM model. The first and (Markevitch 1998> Arnaud & Evrard 1999)j the IC gas 
most noticeable piece of evidence is provided by the fact magg M ^ Mohr ^ aL 19gg) QJ . the gag magg frac _ 
that the X-ray luminosity (L x ) - temperature (T) re- tion f = M /M scale with the IC gas tem per a ture 
lation is steeper than expected from a standard ACDM T diffcren% than expected from pure grav it a tional self- 
model (see, e.g., Arnaud & Evrard 1999, Mushotzky 2003, gimilar evolution of these systems (see Arnaud 20 05 for a 
Arnaud 2005 for a review). This evidence indicated the reyiew) ft hag been suggested that the departures f the 
need for a more detailed look at the IC gas properties. As IC gag gcaling 1&wg from thfi standard se lf- S imil a r model 
a result, the idea that non-gravitational processes can play are Hkely dufi tQ & non _ standard scaling of the mean gas 

indeed a role in determining the cluster structure and evo- , , , A , , u , e ^ ^t 

° density witn temperature and not to a break of tne sen- 

lution has been widely accepted (Evrard & Henry 1991, similarity (see Arnaud 200 5 for a general discussion). The 

Ponman et al. 2003, see Arnaud 2005 for a review). The , , 2/3 , , , , 

, , . ii,! i-i ... . i-i.i i gas entropy, denned as b — KB-in/Pa an d related to the 
picture obtained by the high-sensitivity, high-throughput , , . . ; ., , , 
„ , , /™ , V „„, T J \ , true thermodynamic entropy via a logarithm and an ad- 
X-ray telescopes (Chandra, XMM-JNewton) provided a ,. . . . . 

r . c , . , , r ,, ditivc constant, is considered m this context as a iunda- 
wealth ot precise mlormation on several aspects ot thc 

, , ,. f 1 , j r i mental quantity describing the IC gas and therefore can 

structure and evolution ol clusters and groups ol galaxies. . , J , , , , . 

mi 1 1 1 , • c • i-i i-j.- i- • greatly help to understand thc thermodynamic history of 

I he local relations of various physical quantities (i.e., tern- ° T j; *, Tr ., , , nnnn „„„^n t 

v the IC gas (e.g., Voit et al. 2002, Voit 2005). In the stan- 

~Z ; ~ ~~ ' T „ „ . ~ dard CDM self-similar scenario, thc entropy should scale 
bend offprint requests to: b. Colafrancesco 
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simply as S cx T at any scaled radius. However, it is known 
that the entropy measured at a radius r = 0.1r 2 oo exceeds 
the value predicted by the pure gravitation-based model 
(Pratt & Arnaud 2003, Pratt et al. 2006, Donahue et al. 
2006, Piffarctti et al. 2005), an effect that is especially 
strong for low-T clusters (Ponman et al. 1999, Ponman et 
al. 2003). Various non-gravitational processes have been 
proposed to explain this entropy excess, such as heating 
(from SNs and/or AGNs) before or after the cluster col- 
lapse, or radiative cooling (see, e.g., Voit 2005 for a re- 
view), but none of these is able to provide a conclusive 
explanation of this problem. 

In this framework, recent XMM data (e.g., Pratt & 
Arnaud 2003, Pratt et al. 2006) show a remarkable self- 
similarity in the entropy radial profiles for clusters with 
fcsT ^ 2 keV, a result which is also consistent with a 
stacking analysis of ROSAT data (Ponman et al. 2003). 
The self-similarity of the radial shape of cluster entropy 
profiles proves to be a strong constraint for cluster mod- 
els. For instance, simple pre-heating models, which predict 
large isentropic cores, should be ruled out (Arnaud 2005). 
Furthermore, the slope of the entropy profiles seems to be 
slightly shallower than predicted by shock-heating mod- 
els (Arnaud 2005). A complex interplay of gravitational, 
cooling and galaxy feedback mechanisms is usually consid- 
ered in the attempt to provide a consistent picture (see, 
e.g., Voit 2005 for a review). It is clear that an unambigu- 
ous description of the structure and evolution of groups 
and clusters can be obtained only by studying in detail 
the internal structure of these systems (see, e.g., Arnaud 
2005). 

In this context, there has been no attempt, so far, to 
incorporate a detailed description of the effect of mag- 
netic fields on the IC gas distribution, irrespective of the 
fact that there are several unambiguous and independent 
pieces of evidence for the presence of magnetic fields at the 
^ pG level, especially in the inner regions of galaxy clus- 
ters and groups (see, e.g., Carilli & Taylor 2002, Govoni 
& Feretti 2005) as well as several numerical (see, e.g., 
Dolag et al. 2001) and theoretical motivations (see, e.g., 
Colafrancesco et al. 2005, see also Giovannini 2004 for a 
review) for its presence and extended spatial distribution. 
In this paper we will explore the role of the magnetic field 
B in determining the structural properties of virializcd 
groups and clusters. We will first show in Sect. 2 the effect 
of the magnetic field in determining the properties of the 
IC gas in hydrostatic equilibrium (HE) and their thermal 
structure from the Magnetic Virial Theorem (MVT). We 
will discuss in Sect. 3 the entropy profiles of magnetized 
clusters. In Sect. 4 we will discuss the cluster basic M — T 
relation and we will show in Sect. 5 that the inclusion of 
the B-field provides a flattening of the S — T relation, by 
acting on the IC gas thermal and density properties. We 
will then apply our considerations to the Lx — T relation 
in Sect. 6 and show that the inclusion of the B-field that 
reproduce the S — T relation also recovers the observed 
steepening of the Lx — T relation, from rich cluster to 
poor groups. A systematic exploration of the role of the 



IC gas boundary conditions on the density, entropy and 
X-ray luminosity scaling with the cluster temperature will 
be provided in Sect. 7. We will finally discuss our results in 
Sect. 8 in the context of the existing experimental and the- 
oretical framework, and we will present our basic conclu- 
sions. The relevant physical quantities are calculated us- 
ing Ho = 71 km s _1 Mpc -1 and a flat, vacuum-dominated 
CDM (fi m = 0.3, O a = 0.7) cosmological model. 

2. The structure of magnetized clusters 

We derive in this Section the relevant quantities that de- 
scribe the cluster structure in the presence of a magnetic 
field. These are the density and the temperature of the IC 
gas, both depending, in general, on the magnetic field. 

2.1. Hydrostatic equilibrium of the I CM in the 
presence of magnetic field 

The gravitational potential of a galaxy cluster is deter- 
mined mainly from the DM distribution which contributes 
more that ~ 90% of its total mass. We assume here, for 
simplicity, that the DM density profile is described by the 
NFW profile (Navarro, Frenk & White 1997) 

Pdm = PsVd m {x) = p s ^ ^ 2 , (1) 

where x = r/r s , r s is the scale radius, p s is a normalization 
factor and ydm{%) describes the radial density profile of 
the DM in terms of the adimensional radius x. The scale 
radius r s is linked to the virial radius r„j r by means of the 
concentration parameter c = r V i r /r s , where 

r mr = [M mr /(47rA c (z) i o c (z)/3] 1/3 , (2) 

and p c {z) is the cosmological closure density at redshift z. 
The quantity A c (z) is the average density of the DM halo 
with total mass M V i r = M{< r V i r ) in units of the critical 
density at redshift z = and takes the value A c ~ 100 in 
the reference cosmological model we consider here (Q TO = 
0.3, f^A = 0.7). We adopt here the predicted value of A c (z) 
from Eke et al. (1996). The concentration parameter c 
scales mildly with the cluster mass as 

J M mr V°- 2 

(see Seljak 2000, Colafrancesco et al. 2006). Given the pre- 
vious DM halo structure, the gas density p g and its tem- 
perature T can be obtained from the hydrostatic equi- 
librium (HE) condition once the equation of state for 
the gas is specified. The three cases usually considered 
are: i) the gas follows the DM, p g cx pd m (see, e.g., Wu 
and Xue 2002); ii) the isothermal case T g = const (see, 
e.g., Komatsu and Seljak 2001); iii) the polytropic case, 
T g cx Pg -1 (see, e.g., Komatsu and Seljak 2001, Zhang 
2004). We concentrate here on the isothermal case and we 
improve it by including the effects of the magnetic field. 
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The condition of hydrostatic equilibrium (HE) in the 
presence of a magnetic field B(r) can be written as 



Appendix for details). Thus, the requirement to recover 
the fixed gas mass at the virial radius yields the relation 



dP a dP E 



dr 



GM{< r) 



dr 



-Pg, 



(4) 



where p g = p g (r, B) is the density of the thermal diffuse 
IC gas and P g = P g {r 7 B) is its pressure, with both quan- 
tities depending in general on the magnetic field B. The 
magnetic field pressure, Pb oc B 2 (r) , adds to the left hand 
term of eq.(j4]). The magnetic pressure determines, hence, 
the final distribution of the IC gas density p g (r,B) since 
it tends to counterbalance part of the gravitational pull 
of the cluster thus prohibiting the gas from a further in- 
falling that results in a less concentrated gas core. 
We assume, in the following, a parametric form of the 
magnetic field radial distribution 



B(r) = B* 



Pg( r ) 



10 4 p g {z = 0) 



(5) 



where p g (r) is the solution of the HE equation (see 
Appendix |A|) . p g (z = 0) is the average cosmological den- 
sity of the gas at z — 0, £?» is measured in pG and a is a 
parameter > 2/3 (see, e.g., Zhang 2004, Carilli & Taylor 
2002, Dolag et al. 2001). We assume in our calculations 
here a value a = 0.9 as indicated by numerical simula- 
tions (e.g., Dolag et al. 2001). 

In the Appendix [A] we have derived the solution of the 
HE equation in the presence of a magnetic field B{r) for 
two cases of the equation of state of the IC gas: i) the 
isothermal case and ii) the polytropic case. 
We discuss here specifically the IC gas distribution in 
the case of a isothermal condition in which the Magnetic 
Virial Theorem (MVT) derived for galaxy clusters by 
Colafrancesco & Giordano (2006a) also holds. We will ad- 
dress the polytropic case and more general cases with a 
radial dependence of the cluster temperature elsewhere 
(Colafrancesco & Giordano 2006b, in preparation). 

Under these assumptions, we derived numerically the 
solution of the HE equation for the IC gas density p g (r, B). 
Fig[T] shows the density profile of the IC gas derived from 
the HE condition and normalized to the cluster gas mass 
evaluated at the virial radius as a function of the radius for 
different values of the magnetic field B* . The effect of the 
magnetic field is shown for various masses of the cluster 
M = 0.5,1 and 5 Ms, where Mg is the mass enclosed 
in a 8/i _1 Mpc radius spherical fluctuation at the epoch 
at which it goes non-linear: it takes the value Mg w 2 • 
in our reference cosmological model. The B- 
field provides a pressure that tends to inflate the IC gas 
so that its final radial profile extends outwards, in the 
case of density normalized to the central IC gas density 
in the absence of the B-field. Since observations and our 
modelling of the cluster structure provide a measure of the 
gas distribution for a given mass (within a fixed radius), 
we study the effect of the B-field for a fixed gas mass, i.e., 
extended, say, up to the virial radius, or more generally to 
the radius r& at which the density contrast is A (see the 



,(r,B)=K(M,B)- Pg (r,B = 0), 



with 



K(M, B) 



dxx 2 y g (x,0) 
J C dxx 2 y g (x,B) 



(6) 



(7) 



Such boundary conditions provide a decrease of the cen- 
tral gas density for increasing values of £?* (see FigUJ. 
The effect of B(r) clearly becomes milder for more mas- 
sive systems due to the large DM and gas mass contained 
there. In fact, the normalization function K(M, B) = 
Pg(r,B)/p g (r,B = 0) steeply rises from its low value at 
low masses and rapidly reaches unity at mass scales that 
increase with increasing values of B* (see Figj2j) . Larger 
values of B», hence, increase the cluster mass range that 
is sensitive to the effect of the magnetic field. 

The ratio of the magnetic to IC gas pressure, Pb/P 9 , 
corresponding to the density profiles shown in FigfT] are 
calculated at three different radii r — 0, 0.1r2oo and 
r2oo and are shown in FigfS] as a function of the cluster 
temperature. The ratio between the magnetic pressure, 
Pb oc 8irB 2 (r), and the IC gas pressure, Pg oc p g T g , scales 
as P E /P g oc B 2 p 2 g a ^/Tg ~ BlK 2a - x /T g according to our 
assumption for B(r) (see eqj5|). For relatively high masses, 
when the function K — p g (r, B) /p g (r,B = 0) — > 1 (see 
Fig|2|), the pressure ratio P E lP g oc B 2 T~ X and its ampli- 
tude scales with the value of B 2 . For low masses the gas 
density is strongly affected by the magnetic field and hence 
the function K strongly decreases with decreasing M so 
that the pressure ratio decreases with decreasing M (or 
T) as Ps/Pg oc K 2a -- 1 T~ 1 At an intermediate mass (tem- 
perature) between the low-T tail (where PsjPg strongly 
decreases) and high-T tail (where PsjPg asymptotically 
decreases as 1/T g ), the pressure ratio attains its maximum 
(following the curvature of the function K). The temper- 
ature (mass) location of the maximum of Pb / P g depends 
of the value of B* reflecting the behaviour of the function 
K as a function of -B*, i.e. the behaviour of the gas den- 
sity profiles shown in FigJT] In our specific model for the 
density and temperature of the IC gas in the presence of 
a magnetic field, the value PsjPg is ~ 25% at most and 
peaks at increasing values of T for increasing amplitude of 
The specific values of the maximum of P E /P g , of its 
temperature (mass) location and of the low-T tail depend 
on the specific choice for the exponent a in eq. ([5]) and are 
not universal. Nonetheless, the general trend of the pres- 
sure ratio shown in Fig[3]is preserved for the specific form 
of the radial dependence of the magnetic field as in eq. (0 
and for values of a > 0.5, as indicated by the available 
data. 



2.2. The magnetic virial theorem for galaxy clusters 

The temperature of the IC gas that is used in the HE 
condition is consistently derived, in the framework of 
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Fig. 1. The IC gas density profile as derived from eq.(A19- 
A20) for three different cluster masses (as labelled in the three 
panels), and for values of the magnetic field £>„ = (solid), 0.5 
(dotted), 1 (short dashes) and 3 (long dashes) /iG. The density 
profiles are normalized to the total gas mass at r vir . 



our isothermal model, from the magnetic virial theo- 
rem (MVT). We recall here the main aspects of the 
IC gas temperature derivation from the MVT following 
Colafrancesco & Giordano (2006a). 

Under the assumption of an ICM in hydrostatic equilib- 
rium with the potential well of a spherically-symmetric, 
isolated, virialized and magnetized cluster, the general re- 
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Fig. 2. The function K(M,B) = p g (r, B)/ p g (r, B = 0), eval- 
uated at r = 0.1r2oo, is shown as a function of the cluster mass 
and for three different values of the magnetic field B, = 0.5 
(cyan), 1 (red) and 3 (green) pG. 



lation between the ICM temperature T g and the cluster 
virial mass M V i r is obtained by applying the MVT that, 
for a static and isothermal galaxy cluster, reads 



2K + 2U + U B + W = 0. 



(8) 



Here Ub is the magnetic energy of the system, U is the ki- 
netic energy of the gas (we assume that the Dark Matter 
is cold and collisionless), K is the Dark-Matter particle 
kinetic energy and W is the total potential energy (see 
Colafrancesco & Giordano 2006a for details). The previ- 
ous eq. ([8|) holds specifically in the absence of an external 
medium. For the more general case of a cluster that is im- 
mersed in a Inter Galactic Medium (IGM) or in an exter- 
nal medium that exerts an external pressure P ex t, eq.© 
yields the formula for the temperature of the gas in virial 
equilibrium 



pm p 



£,G M v 



1 



Ml 
Ml 



P, 



p 

1 V. 



(9) 



where P V i r = 
given by 



-) 1 with usually £ ~ 1, and M$ is 



M$ ~ 1.32 • 10 13 M, 



1(c) 



1/2 



JtGJ \MpcJ 



(10) 



0)) 2a x 2 y 2a (x,B 



Here 1(c) = ^(p g (r = 0)/p g (z = U)f-x"y' g 
0)dx, and y g (x,B = 0) = p g (x)/p g (x = 0) is the gas 
density profile normalized to the central gas density (i.e. 
the solution of the hydrostatic equilibrium equation in the 
absence of a magnetic field, see Appendix for details) . 
For the case P ex t = and B — 0, one finds = and 
the well-known relation 



k B T g (B = 0) 



^pbrripW 



(11) 



S. Colafrancesco & F. Giordano: Structure and evolution of magnetized clusters 




For B > 0, one finds > and the temperature of the 
gas at fixed M vir reads 
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Fig. 3. We show the ratio Pb/P 3 between the IC gas pressure 
and the magnetic pressure as a function of the cluster temper- 
ature for different values of the magnetic field £?* = 1 (red), 3 
(green), 5 (blue) and 10 (cyan) /xG. The ratio Pb/-P 3 is shown 
at three different cluster radii: r = (upper panel), r = 0.1r2oo 
(mid panel) and r — r2oo (lower panel). 



kT g = kT g {B = 0) 1 - 



Ml 



Pe, 



(12) 



The temperature of a cluster calculated in the presence of 
a magnetic field (with P ext = 0) is lower than that given 
by eq. fPPj) because the additional magnetic field energy 
term U B adds to the MVT. 

The presence of an external pressure P ex t tends to com- 
pensate the decrease of T g induced by the magnetic field. 
In the following, we consider values of P ex t (which are rel- 
evant for the structures we consider here) as estimated 
in the IGM and in the rich cluster outskirts. For values 
of the temperature, Tjgm ~ 10 6 K, and density, tligm ~ 
10 _5 cto -3 , as estimated by the WHIM structure around 
large-scale overdensities (see, e.g., Fang & Bryan 2001), 
P ext ~ 1.7 • 10" 3 eV cm- 3 (n /GM /10- 5 cm- 3 )(T /GM /2 • 
lO 6 !^). However, in the outer regions of massive clusters 
(at r £ r V i r , where the main source of pressure could 
be due to the momentum of the infalling gas accreting 
onto the cluster), studies of the mean projected temper- 
ature profile for the cluster sample derived by Piffaretti 
et al. (2005, see their Fig. 4) indicate that the exter- 
nal gas pressure can reach values up to P ext ~ 0.2 eV 
cm~ 3 (n/KT 4 C7n~ 3 )(T 9 /1.7 • 10 7 K) (we consider here a 
typical rich cluster with T g — 10 keV). In the case of such 
a hot cluster, the value of P ex t is a significant fraction, 
~ 4%, of the central ICM pressure and ~ 50% of the ICM 
pressure at the virial radius for a typical cluster. However, 
for lower-temperature clusters, such values of P ex t are ex- 
pected to be sensitively lower. A value P ex t ~ 0.2 eV 
cm~ 3 , as estimated at the outskirts (r Si r V i r ) of rich, 
hot clusters, can be considered, hence, as an upper bound 
to Pext] hi addition, an exact determination of the total 
cluster mass (which is subject to various systematic uncer- 
tainties, see, e.g., Rasia et al. 2006) requires an extension 
beyond r V i r . Lower values of P ex t, down to the value found 
in the WHIM, have progressively less importance for the 
MVT. 

For reasonable values of ^ a few /LtG, the quantity 
M? > M% ir ■ (Pext/Pvir) in eq.© and the main effect 
is a reduction of the cluster temperature which is more 
pronounced for less massive systems, where becomes 
comparable to M V i r . The effects of the magnetic field and 
of the external pressure are greater for low-M clusters (see 
Colafrancesco & Giordano 2006a and Figj5j) . 



re-obtains (here fj, — 0.63 is the mean molecular weight, 
corresponding to a hydrogen mass fraction of 0.69, m p is 
the proton mass and k B is the Boltzmann constant). 



3. The radial entropy profiles of magnetized 
clusters 

Given the previous results on the temperature and on the 
gas density in the presence of a B-field, we can now discuss 
the impact of the magnetic field on the entropy profile of 
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galaxy clusters. 

The entropy of the IC gas is usually defined as 
q _ k sT g 

J 2/3 

Pg 

and, in our model, it depends on the magnetic field 
through the B-dependence of the cluster temperature 
T g (B), as derived from the MVT (see eqE]), and from the 
IC gas density p g (r, B), as derived from the HE condition 
(see eqs.Al9-A20). 

Our model provides a simple but consistent description of 
the B-dcpcndcnce of the (isothermal) cluster entropy and 
contains all the relevant physical information. However, 
the detailed effect of the magnetic field on the cluster 
entropy should be addressed in a truly magneto-hydro- 
dynamical (MHD) model of cluster formation and evolu- 
tion, which is able to catch the details of the turbulent 
amplification and of the evolution of the B-field associ- 
ated with the shocked gas in the cluster environment. 
Such a task is well beyond the scope of this paper and 
could be better tackled through MHD cosmological sim- 
ulations and/or semi- analytical methods that will be ad- 
dressed elsewhere. Nonetheless, the model we present here 
is able to describe the overall properties of cluster entropy, 
namely its radial behaviour and its correlation with the 
cluster temperature. In Fig|4]we show the radial entropy 
profiles of the IC gas calculated by using the IC gas den- 
sity profiles normalized to the cluster mass gas (these are 
shown in FigH]). The entropy profiles we obtain flatten in 
the inner cluster region for increasing values of the mag- 
netic field . This behaviour reflects the decrease of the 
central IC gas density for increasing values of the magnetic 
field, as shown in FigQ] The bottom panel in Fig[4] shows 
also the comparison of our entropy profiles for a typical 
M V i r = 5M§ cluster with the region enclosed by the mean 
plus/minus the 1 er standard deviation of the scaled en- 
tropy profiles for ten clusters observed with XMM (Pratt 
et al. 2006). These clusters have, in fact, quite high masses 
around the reference value we consider here. The predic- 
tions of our model agree quite well with the data, and the 
widening of the observed entropy profiles at small radii 
could reflect the dispersion of the B-dependent density 
profiles found for high values of B r . In this respect, the 
XMM cluster data, taken at face value, limit the ampli- 
tude of the magnetic field to values £?* ^ 7 fxG for this 
mass range. 

The mean temperature used to investigate the scaling 
properties of the XMM clusters has been estimated in the 
region 0.1r 2 oo < T < 0.5r 2 oo and it is quite constant for 
most of the XMM clusters (see Pointecouteau et al. 2005). 
This fact allows us to directly compare the predictions of 
our isothermal model with the observed entropy profiles 
derived by Pratt et al. (2006). 

Furthermore, our model predicts a sensitive flattening of 
S(r) at smaller radii < 0.2r2oo which is more evident for 
higher values of This behaviour is similar to the en- 
tropy flattening shown in the inner cool cores of the clus- 
ters studied with Chandra (see, e.g., Donahue et al. 2006). 
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Fig. 4. IC gas entropy profile for various values of the mag- 
netic field B t — 0, 0.5, 1 and 3 ^iG and for three different cluster 
masses: 0.5Mg (top panel), Ms (mid panel) and 5M% (bottom 
panel). The density profiles are normalized to the total gas 
mass at r vir . The shaded region in the last plot for a cluster 
with Mvir = 5Ms show the 1 sigma area spanned by the nor- 
malized entropy profiles derived by Pratt et al. (2006) for 10 
nearby clusters. 



Note, however, that a closer comparison of the (isother- 
mal) model we consider here with these cool cluster cores 
(where one observes temperature jumps of the order of 
a factor of 2-3 with respect to the outer temperature, 
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see, e.g., Colafrancesco, Dar & DeRujula 2004 and ref- 
erences therein) cannot be done straightforwardly, and it 
demands a more detailed modeling of the magnetized core 
regions that will be presented elsewhere (Colafrancesco & 
Giordano 2006b). 

4. The M-T relation for magnetized clusters 

To derive consistent structural relationships for clusters 
and groups, we first need to derive their M — T rela- 
tion. The M — T relation is a basic structural relation for 
groups and clusters and it is the starting point to study 
the entropy and X-ray luminosity scaling with T and/or 
M. The T — M relation for magnetized clusters has been 
derived by Colafrancesco & Giordano (2006a) and is also 
shown in FigJS] where it is compared to the data for the 
XMM clusters used to derive the entropy profiles previ- 
ously discussed (see Pratt et al. 2003). Our isothermal 
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Fig. 5. The T 3pec t r — M200 relation at z = for clusters that 
contain a magnetic field B t — (black), 10 (magenta) and 
30 fiG (red), and for values of P ex t = (dashed curves) and 
Pext = 0.2 eV cm -3 (solid curves). Data are taken from Arnaud 
(2005). 



model is consistent with the previous data because these 
are obtained by using an overall spectroscopic tempera- 
ture (which is constant in the radial range (0.1 — 0.5)r2oo 
and also because the values of T spec t r , have been obtained 
by performing an isothermal fit of the cluster spectra ob- 
tained in the same (0.1 — 0.5)r 2 oo radial range. 
The normalization of the observed T spectr . — M200 rela- 
tion is however found to be discrepant by ~ 30% with 
respect to the value derived from numerical simulations 
including only gravitational heating (see, e.g. Evrard et 
al. 1996), a well known problem (see, e.g., Arnaud 2005 
for a review) which usually requires one to re-normalize 



the predictions of analytical models. In our model we nor- 
malize the T spectr . — M200 relation for the case B = to 
the observed data derived by Arnaud et al. (2005) by as- 
suming T spectr = T g , where T g is obtained from eq.©, 
and a temperature boost of a factor £ as 1.5 in eq.©, 
as can be expected from the continuous shock-heating of 
the IC gas within the virial radius after the formation of 
the original structure (see, e.g., Makino et al. 1998, Fujita 
et al. 2003, Ryu et al. 2003). Such a prescription for £ 
has been used by the previous authors for unmagnetized 
clusters in the absence of external pressure. A systematic 
effect which tends to increase the cluster temperatures is 
the presence of a minimal external pressure in eqs.(9) and 
(12). A value P ext ~ (0.1 - 0.2)P mr (like that found in 
the IGM around clusters) could easily accommodate for 
an overall boost value of £ » 1.5 x {P ex t/P V ir) > 1-5 and 
reasonably in the range 1.65 — 1.8, for the previous val- 
ues of P ex t ■ Given the large theoretical uncertainty on the 
non-gravitational heating efficiency, we adopted an over- 
all value £ ~ 1.8 to normalize our prediction to the data 
point at Af2oo ~ 3-Ms in Fig. 5, which is the point with the 
smaller intrinsic error. However, values of £ in the plausi- 
ble range 1.5 — 1.8 marginally change our predictions. 
The relation M200 — -77 M V i r is also found in our mass 
scale definition. 

Small variations of temperatures with respect to their un- 
magnetized values are found for massive clusters since 
M V i r in this mass range and the value of P ex t 
has little or negligible effect (see Fig[5]). This is in agree- 
ment with the results of numerical simulations (Dolag 
et al. 2001a). However, when becomes comparable 
to M V i r , the IC gas temperature becomes lower than its 
unmagnetized value and the T — M relation steepens in 
the range of less massive systems like groups and poor 
clusters. The temperature T g formally tends to zero when 
M,p — > M V i r (l + Pext/Pvir) 1 ^ 2 ■ However, this limit is un- 
physical since it corresponds to an unstable system in 
which the magnetic pressure overcomes the gravitational 
pull. Thus, any physical configuration of magnetized viri- 
alized structures must have < M vir (l + P ext / Pvir) 1 ^ 2 ■ 
The effect of P ex t counterbalances the effect of the S-ficld 
on the T — M relation, and its amplitude increases for low- 
M systems, (see Figj5] see also Colafrancesco & Giordano 
2006a for details). 

5. The S — T relation for magnetized clusters 

Using the T — M relation presented in Sect[4l we can now 
derive our predictions for the S — T relation for clusters 
with a magnetic field. Here, we will compare our models 
for the IC gas distribution in the presence of a magnetic 
field with the data of the S — T relation (see Ponman et 
al. 2003, see also Arnaud 2005 for a review). The entropy 
in this correlation has been evaluated at a distance 0.1r2oo 
from the cluster center and thus we have to evaluate both 
the IC gas density and its temperature at this radius. 

Due to the decrease of the IC gas density for increasing 
values of the magnetic field (see Fig [J), the density of the 
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IC gas evaluated at a distance 0.1r 2 oo exhibits a peculiar 
trend with increasing cluster temperature (mass) as shown 
in Fig[5J The presence of a temperature dependence of the 





k B T Ilr (keV) 



Wn» (keV) 

Fig. 7. The entropy S(0.1r2oo) as a function of the temper- 
ature for various values of the magnetic field B* as labelled. 
The binned entropy data are taken from Ponman et al. (2003). 



Fig. 6. The temperature dependence of the IC gas density 
evaluated at 0.1r2oo for various values of the magnetic field 
B* = (black curve), 0.5 fiG (red), 1 fiG (yellow), 2 ^G 
(green), 3 /iG (blue), 5 fiG (magenta) and 10 fiG (cyan). 



gas density p g (0.1r 2 oo), even in the case B = 0, is a con- 
sequence of the HE condition (see eq.4) , independently of 
the chosen equation of state for the gas in the MVT. The 
specific solution of the HE condition for the unmagnetized 
case (B = 0) is given in eq.(A.lO) of the Appendix. It is 
clear that p g (r) oc exp[— 3c/m(c) j Q dum{u)/u 2 ] exhibits 
a temperature dependence (like that shown in Figj6]) since 
the concentration parameter c depends on T g through 
eq.(3). 

For each fixed value of B*, the presence of a B-field pro- 
duces a sharp rise of the gas density p g (0.1r 2 oo) in the 
low-T region followed by a smooth decline at higher tem- 
peratures where the density tends to the unmagnetized 
value (solid black curve). This effect is produced by the 
physical boundary condition for the IC gas density, i.e. by 
the requirement to normalize the IC gas density to the 
total gas mass at r vir for each cluster. The sharp rise of 
p s (0.1r 2 oo) with the cluster temperature (mass) is pro- 
duced by the effect of the function K(M,B) (see FigEJ: 
when K — > 1 at high masses, then p g (0.1r 2 oo, B) tends to 
its value for unmagnetized clusters with B = (i.e., the 
black curve in FigJSl see also discussion in Sect. 2.1). 
Given the temperature behaviour of p 9 (0.1r 2 oo) and the 
value of the IC gas temperature at the same radius, the 
S — T relation for magnetized clusters predicted by our 
model can be calculated using eq. (ll3p and is shown in 
FigO We stress that, based on the HE condition for the 
IC gas, one has to expect an entropy scaling with the gas 



temperature which is different from the one, S ~ T, usu- 
ally associated with "standard self-similar" models, even 
in the case B = 0; the actual scaling is, in particular, 
steeper than S ~ T due to the decreasing behaviour of 
p s (0.1r 2 oo) with increasing T (see FiglBJ). 
The effect of a fixed amplitude of the magnetic field -B* 
(set constant here for any cluster) is to flatten the S — T 
relation due to the combination of the decrease in T g 
caused by the MVT and to the decrease of the IC gas 
density p 9 (0.1r 2 oo) caused by the HE condition. The en- 
tropy inversion shown for low-T (at fixed value of -B*) in 
Fig|7] reflects the sharp cutoff at low-T in the gas den- 
sity p g (0.1r 2 oo) shown in Fig[5J The binned entropy data 
derived from Ponman et al. (2003) shown in Fig[7]are en- 
closed by the two curves for B* = 2.5 and 7 fiG. This 
indicates that a magnetic field whose amplitude decreases 
with decreasing cluster temperature (mass) as £?» ~ T^ q 
with rj « 0.5 may reproduce, on average, the whole dis- 
tribution of the data. The best-fit trend of the entropy 
S(T) evaluated by assuming this phenomenological scal- 
ing of the magnetic field amplitude, £>» ~ T 5 , is shown 
in FiglU We also report in Figl9] the best-fit relation for 
the amplitude of the magnetic field, B* oc T 20 q, plus its 
variation allowed by the IC gas entropy data. 

We notice that the effect of a magnetized IC gas pro- 
duces an inversion of the entropy at low temperature val- 
ues, i.e. in the group region, as the one shown by the 
unbinned data derived from Ponman et al. (2003). FigfTUl 
shows our best-fit prediction plus the uncertainty region 
for the S — T relation compared to the unbinned entropy 
data of Ponman et al. (2003, see their Fig. 4). We also plot 
in Fig [11] the S — M relation predicted by our model. 
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Fig. 8. The entropy S^O.l^oo) i s shown as a function of the Fig. 10. Same as Fig(S]but for our best-fit model compared 



temperature for a magnetic field scaling B t ~ T^ 00 with r\ = 0.5 
(solid curve). The dashed area enclosed by the two entropy 
curves calculated for B„ ~ T%qq and for _B, ~ T^gg shows the 
range of variation for B, allowed by the binned entropy data. 
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to the unbinned data distribution of Ponman et al. (2003). 




M /M. 



1 10 

k B T z „ (keV) 

Fig. 9. The function B* = 2.6/iG(T 20 o /keV) 05 which best fits 
the S — T relation in Fig[8] is shown by the solid curve. The 
shaded area shows the range of variation of B*(T2oo) enclosed 
by the curves B, = 2AfiG(T 20 o/keV) OA (lower bound) and 
= 2.9nG(T2oo/keV) (i (upper bound). 



6. The Lx — T relation for magnetized clusters 

We show in this section that the same physical model that 
is able to recover the flattening of the S — T relation also 
provides an explanation for the steepening of the Lx — T 
relation, from the scale of clusters down to the group scale. 



Fig. 11. The entropy S^O.lraoo) as a function of the cluster 
virial mass M V ir (in units of Ms) for various values of the 
magnetic field B, as in Fig[7] 



Since the X-ray emissivity is mainly provided by thermal 
bremsstrahlung emission 

e brem w 3 • 10- a7 ergs~W 8 p2(r, B)T^/ 2 (B) , (14) 

(here we adopt its expression for ksT g ^2 keV, see 
Sarazin 1988), the bolometric X-ray luminosity of a cluster 
with mass M is given by 

L X o,T^ 2 (B)pl(0,B)C^\ dxx 2 y g (x,B) (15) 
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which reads, in practical units, as 

' k B T g (B) 
AkeV 



L x w 2.8 • 10 45 ergs 



1/2 



10 3 cm 3 



1 (M„ 



Ms 



dxx 2 y g (x, B) . 



(16) 



The presence of a magnetic field provides a strong decrease 
of the cluster Lx at low temperature due to the decrease of 
the central IC gas density (see Sect. 2) and also a decrease 
of the cluster temperature from the MVT with respect to 
their unmagnetized values. 

The result of a simple calculation based on the HE, on the 
MVT and on the expression of Lx as given in eqs. (|15I16[) is 
reported in FigfT2land it shows that the effect of the mag- 
netic field is to steepen the Lx~T relation. We stress that 




k„T (keV) 

Fig. 12. The bolometric X-ray luminosity Lx as a function 
of the gas temperature for various values of the magnetic field 
B* — (black curve), 2.5 /xG (yellow), 3.7 fiG (green), 5 jj,G 
(blue), 6 fj,G (red) and 7 /iG (cyan). Data are taken from 
Arnaud & Evrard (1999, blue dots), Allen et al. (2001, red 
dots) and Ponman et al. (2003, magenta dots). 



the X-ray luminosity data shown in Fig |12l are enclosed by 
the same curves that enclose the S — T data in the same 
temperature range. Moreover, as in the S — T relation, 
a magnetic field whose amplitude decreases with decreas- 
ing cluster mass as £?* ~ j^-5±o.i j s a j-,i e to re p rocm ce 
the whole distribution of the Lx — T data from groups to 
clusters. The best-fit trend of the Lx — T relation eval- 
uated by assuming this phenomenological scaling of the 
magnetic field amplitude which fits the entropy data (see 
Figs. 8 and 10) is shown in Fig[T3J 

We also show in FigfTH the Lx — M relation that is 
consistent with the Lx — T relation previously discussed. 

The fact that the same model for magnetized clusters 
which is able to reproduce the S — T relation (see Sect. 5) 




k B T (keV) 

Fig. 13. The bolometric X-ray luminosity Lx as a function 
of the gas temperature for a magnetic field scaling _B* ~ T^ 00 
with r\ = 0.5 ± 0.1. Note that this magnetic field scaling is the 
same that fits the S — T relation, as shown in Fig|8] 




M Ilr /M a 



Fig. 14. The bolometric X-ray luminosity Lx as a function of 
the cluster mass for various values of the magnetic field B* = 
(black curve), 2.5 fj,G (yellow), 3.7 piG (green), 5 [iG (blue), 6 
/xG (red) and 7 /iG (cyan). 



also recovers the Lx — T relation is not a coincidence, 
but is a consequence of the (reliable) physical descrip- 
tion of the X-ray cluster structure. In fact, the entropy 
S = T g l p 2 J^ calculated at a scale 0.1r 2 oo (he., within the 
cluster core) sets the values of the gas density and tem- 
perature at the same spatial scale, and hence the value of 
the X-ray luminosity Lx oc p 2 g T g , which is mostly deter- 
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mined by the cluster core properties. Therefore, adjusting 
the parameters of our model to match the observed S — T 
relation implicitly ensures that the resulting Lx — T rela- 
tion is close to the observed one. 

7. Exploring the effects of boundary conditions 

In the previous sections we have shown the effect of the 
magnetic field on the S — T and Lx — T relations for the 
case in which i) the external pressure was negligible (the 
formal case P ex t — has been adopted in Sects. 4-6) and 
ii) the boundary condition used to normalize the IC gas 
density profile given by eq.©, i.e. a constant gas mass 
within the cluster virial radius. 

For completeness, we describe in this section the effects 
on cluster temperature, entropy and X-ray luminosity in- 
duced by changing the previous boundary conditions. 

7.1. Normalizing to the cluster virial pressure 

Normalizing the IC gas density profile to have a constant 
pressure at the virial radius, P g (r V i r , B) = const., yields a 
decrease of the central gas density with increasing values 
of Btf which is larger than when normalizing the same den- 
sity profile to M g (r v i r ); this happens because, for a given 
cluster mass (temperature) , the former condition requires 
Pg{ r vin B) = p g (r viri B = 0). This fact implies a be- 
haviour of the function K(M, B) = p g (r, B)/p g (r, B = 0) 
that is slightly different from that shown in Figf21 and con- 
sequently a slightly different behaviour of the scaling of 
/ o g (0.1r2oo) with the cluster temperature (mass) as shown 
m Fig [H Since the maxima of the function jO g (0.1r2oo) 
are moved to lower temperatures in this case, the inver- 
sion of the S — T relation is also moved at lower temper- 
atures and the flattening of the Lx — T relation is found 
at lower values of the X-ray luminosity and lower tem- 
peratures (see Fig fTHl) . These results show that both the 
entropy inversion and the flattening of the Lx—T relation 
at low temperatures are a realistic effect to be expected in 
models of magnetized clusters; moreover, these results do 
not strongly depend on the adopted boundary conditions 
for the gas density profile of the clusters. This reinforces 
our conclusions on the ability of models of magnetized 
clusters to reproduce the basic structural properties of X- 
ray clusters. 

To conclude our exploration of the effects of boundary 
conditions, we will address in the next section the effects 
on gas density, entropy and X-ray luminosity caused by 
changing the external pressure P ex t at the cluster bound- 
ary. 

7.2. Changing the cluster external pressure 

Changing the value of the external pressure P ex t in the 
MVT modifies the cluster temperature and also alters the 
solution p g (r, B) of the HE condition as shown in FigfTTl 
In fact, increasing the value of P ex t yields an increase of 



the overall cluster temperature (see eq|9j) and hence a de- 
crease of the central IC gas density, to be consistent with 
its boundary conditions (see FigfT7|. The increase of P ex t 
in the MVT and in the HE condition causes a larger de- 
crease of the IC gas density in the cluster core (even in 
the absence of a magnetic field) and hence a further bend- 
ing of the temperature dependence of /9 9 (0.1r2oo) in the 
low-T region than the case in which P ext = (see FiglBJ). 
This behaviour implies that the S — T relation is flatter 
at low-T with respect to the case in which P ext — 0, even 
for £>„ = (see Fig|T8l). In addition, the entropy inver- 
sion which is produced in the low-T region by the effect of 
magnetic field is further enhanced by the increasing value 
of P ex t- As a direct consequence of this fact, the Lx —T re- 
lation computed for values P ex t > becomes steeper than 
that derived in the case P ex t = 0, even in the absence of a 
magnetic field. The presence of a magnetic field (like the 
one B» oc T 5 that best fits the S — T and the L x - T 
relations for P ext = 0, see Figs 151 and [T3|) further increases 
the steepness of the predicted Lx — T relation towards 
low temperatures, thus requiring lower values of the mag- 
netic field to reproduce, in this case, the distribution of 
the available data (see Figll8[). 

Finally, the combination of large values of P ex t and a 
different boundary condition (like the one requiring us to 
normalize the IC gas density profile to the virial pressure, 
that we previously explored in this section) has the com- 
bined effect of further increasing both the flatness and the 
inversion of the S — T relation and, consistently, the steep- 
ness of the Lx — T relation in the low-temperature region 
(see FigfT^j) . We stress that in this last case the S — T 
and the Lx — T relations evaluated for a null magnetic 
field are no longer representable by simple power-law be- 
haviours, since the presence of a positive external pressure 
alters the cluster temperature derived from the MVT and 
the value of p fl (0.1r2oo)- This suggests that, in this case, 
the observed Lx — T relation (as well as the S-T relation) 
could be roughly reproduced by the combination of high 
values of P ex t with a more modest contribution from the 
magnetic field (with respect to the case P ext = 0), because 
the effect of P ext in the MVT and in the HE conditions of 
our model (see Sect. 2) tends to mimic - at a given bound- 
ary condition for the IC gas at the cluster virial radius 
- that of an increasing magnetic field on the IC gas den- 
sity profile. In fact, the external pressure compresses the 
IC gas and tends to rise its temperature, thus requiring 
a decrease of the core density in order to fulfill the IC 
gas boundary condition at r v i r . Both the boundary condi- 
tions that have been considered here (M g (r V i r ) = const., 
or Pg(r V i r ) — const.) provide hence a decrease of the IC 
gas density in the cluster core which affects both the en- 
tropy and the X-ray luminosity, as we have shown here. 
By stretching our analysis to the extreme conclusion, we 
could say that the S — T and the Lx — T relations could 
be reproduced with an appropriate distribution of exter- 
nal pressure values, decreasing from the rich cluster scale 
down to the group scale, with only a small effect due to 
the unavoidable presence of a low-amplitude intra-cluster 
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Fig. 15. Left. The radial density profiles of the IC gas normalized to the same virial pressure P V i r for values of the magnetic 
field B* — (solid) and B* = 3 fj,G (dashed cyan) are compared to the analogous density profiles normalized to the total gas 
mass at rw (dashed black). A cluster with M — Ms is considered here. Right. The temperature dependence of the IC gas 
density evaluated at 0.1r2oo for the same values of the B-field as in Fig|6]but for a normalization of the IC gas density to the 
virial pressure P V i r - 




Fig. 16. Left. The S — T relation for clusters with IC gas density profiles normalized to the virial pressure P V i r is shown for the 
same values of the magnetic field B* as in Fig[7] Right. The Lx — T relation for clusters with IC gas density profiles normalized 
to the virial pressure P V i r is shown for the same values of the magnetic field B, as in Fig ll2l 



magnetic field. However, we recall that values oiP ext « 0.2 
eV cm -3 (considered in FigsfT7l-fT9|) are found at the out- 
skirts (r r v i r ) of rich, hot clusters (with mean temper- 
atures <; 10 keV, see Sect. 2. 2) and they are, therefore, 
extreme values of the external pressure which are much 



higher than those usually found around low-T clusters and 
groups. This means that the temperature scaling of the 
cluster entropy and X-ray luminosity (with and without 
a magnetic field) for such a high value of P ex t considered 
as representative of the whole population of clusters (from 
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Fig. 17. Left. The radial density profiles of the IC gas evaluated in the presence of an external pressure P ex t = 0.2 eV cm -3 
(solid curves) are shown for values of _B« =0 (black) and B, — 3 fiG (green) and are compared to the analogous profiles 
evaluated with P ext = (dashed curves for the same magnetic field values). A cluster with M = Ms has been chosen. Right. 
The temperature dependence of the IC gas density evaluated at 0.1r2oo is shown for the same values of the B-field as in Fig[6] 
but for a value of the external pressure P ext = 0.2 eV cm -3 . 





Vax, (keV) 
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Fig. 18. Left. Comparison of the best fit S — T relation evaluated for P ex t = (dashed cyan, see Fig(8]) and the one calculated 
for the same parameters but for P ext = 0.2 eV cm -3 (solid cyan). The S — T relations evaluated in the case B* = with 
Pext = (dashed black) and B, = with P ext = 0.2 eV cm -3 (solid black) are also shown for comparison. Right. The same 
cases shown for the S — T relation in the left panel of this figure are shown here for the Lx — T relation. Different curves refer 
to the same cases as discussed in the left panel. 



rich ones down to small groups) has to be considered un- 
likely; a lower value of P ex t closer to the value found in 
the IGM is better suited to the case of low-T systems 



(see, e.g., our discussion in Sect. 2, see also the results of 
the temperature and density profile analysis of a sample 
of poor and rich clusters studied by Piffaretti et al. 2005), 
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Fig. 19. Left. The S — T relation for clusters with IC gas density profiles normalized to the virial pressure P„j r and with 
Pext = 0.2 eV cm -3 . Curves are for the same values of the magnetic field B* as in Fig[7] Right. The Lx —T relation for clusters 
with IC gas density profiles normalized to the virial pressure P v i r and with P ext = 0.2 eV cm -3 . Curves are for the same values 
of the magnetic field B, as in Fig |12l 



and thus the results presented in Sects. 3-6 seem to be 
more reliable. 

The results discussed in this Section have been pre- 
sented for illustrative purposes and are based on a para- 
metric variation of the leading physical quantities and 
boundary conditions. A more detailed exploration of these 
effects requires a thorough analysis of the density and tem- 
perature profiles of a large sample of groups and rich clus- 
ters out to radii r <; r V i r , an analysis which is crucial for 
this issue, but it is not yet available. The previous anal- 
ysis should be combined with the complementary study 
of the evolution of the intra-cluster magnetic field with 
the cluster scale (mass, radius) which has just begun and 
that will receive a substantial boost from the advent of the 
next generation radio observatories with high sensitivity 
and spatial resolution (e.g, SKA, LOFAR). Given the lack 
of detailed data on these subjects, we will tackle in detail 
the theoretical aspects related to these issues in a forth- 
coming paper (Colafrancesco et al. 2006, in preparation). 

8. Discussion and conclusions 

Our results provide an analytical description of the effects 
of a magnetic field on the density structure and on the en- 
tropy profiles of galaxy clusters in the context of the HE 
condition and of the MVT. We have shown that such a 
description is able to provide, at once, a possible explana- 
tion of three problematic aspects of the cluster structure: 
i) the flattening of the entropy profile in the cluster cen- 
ter (together with the possible entropy inversion near the 
center of the structure); ii) the flatness of the S — T re- 



lation (together with the entropy inversion at the scale of 
the galaxy groups); iii) the increasing steepening of the 
Lx — T relation from the cluster scale towards the group 
scales. 

The inclusion of the effects of the magnetic field in the 
description of the intra-cluster gas structure is motivated 
by the increasing evidence for the presence of a magnetic 
field in these cosmic structures (see, e.g., Carilli and Taylor 
2002, Govoni and Feretti 2005 for recent reviews). 

In the framework of our model, all the available data 
indicate that an increase of the magnetic field ~ 
T°' 5±01 is able to reproduce, at the same time, both the 
S — T and the Lx — T relations. Such a scaling is a rea- 
sonable expectation of hierarchical clustering scenarios in 
which the magnetic field is coupled to the amount of gas 
collected in the deepening potential wells of increasing 
mass systems. The available evidence on the amplitudes of 
cluster magnetic fields for individual clusters and groups 
seem to be consistent with this result (see also Cassano et 
al. 2006 for an independent study). 

Since the theoretical results presented here partly depend 
on the assumed boundary conditions for the IC gas at 
r « r v ir, we have also provided a systematic exploration 
of the effects of changing the IC gas boundary conditions 
on the structure and evolution of X-ray clusters. Such an 
analysis has been able to show that the study of the mag- 
netic field effects on the IC gas must be addressed in full 
complementarity with the pressure and density study of 
the cluster external regions at r ^ ?w. 

How does the model and the results we have presented 
here compare to other theoretical explorations of the IC 
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gas physics? 

It has been suggested (see, e.g., Arnaud 2005 for a review) 
that the steepening of the Lx — T relation could be due 
to a systematic change (increase) of the mean gas density 
with T (i.e., a simple modification of the scaling laws) or 
to a variation of the cluster shape with T (i.e., a breaking 
of self-similarity) . Many studies have attributed the devia- 
tions from self-similarity to an early episode of preheating 
by supernovae or active galaxies so that the heat input 
that cstabilishes a uniform minimum entropy level in the 
intergalactic medium breaks the self-similarity because the 
extra entropy makes the gas harder to compress as it ac- 
cretes into dark matter halos (e.g., Kaiser 1991, Evrard 
& Henry 1991). However, several analyses suggested that 
the energy input required to explain the observed relations 
through global preheating seems to be quite extreme (see, 
e.g., Donahue et al. 2006, Arnaud 2005). Alternative ex- 
planations (Voit & Bryan 2001) have argued that the en- 
tropy scale responsible for the breaking of self-similarity 
is not a global property of the IGM but rather a condition 
set by the radiative cooling. In such a scenario, cooling and 
feedback seem to conspire to deplete the amount of gas be- 
low the cooling threshold: thus, some of the low-entropy 
gas condenses and the feedback subsequently raises the 
entropy of the remaining gas until both cooling and feed- 
back shut down. Variations on this theme (see, e.g., Voit 
et al. 2002) are also able to provide models that may ac- 
count for several X-ray observables (like the Lx — T and 
the T — M relations, the X-ray surface brightness pro- 
files of clusters and their temperature gradient). However, 
these models are not yet able to explain why clusters differ 
in the amount of gas that still remains below the cooling 
threshold. 

In such a complex theoretical framework, some recent 
studies of the self-similarity in the cluster gas tempera- 
ture and entropy profiles suggest that the departures of 
the entropy and X-ray luminosity scaling laws from the 
standard self-similar model could be due to a change in 
the temperature (mass) dependence of the normalization 
of the gas density profiles (see Arnaud 2005 for a general 
discussion). 

In our model a non-standard scaling of the IC gas den- 
sity with the cluster temperature (mass) is provided (see 
FigsQ] and [6]) by the inclusion of the magnetic field ef- 
fects on the HE (and the MVT) condition of the IC gas, 
which then changes both the radial distribution of the IC 
gas (for a given cluster mass and B-field value) and the 
relation between the gas temperature and cluster mass. 
It follows that a consistent description of the magnetized 
ICM can provide a simple explanation of several (or of all) 
of these still open problems and thus weakens the need for 
the inclusion of non-gravitational effects - like e.g., shock 
heating, AGN or SN feedback - which have been proposed, 
so far, for the explanation of these features. 

Our study (far from being conclusive) shows that a 
more global description of the IC gas physics may help in 
solving some of the crucial aspects of cluster structure and 
evolution. 



In these respects, we want to stress, in the following, some 
of the limits of this initial study. 

The isothermal case that we explored here in detail is 
not valid for all clusters and groups, even though it seems 
to be a valid approximation in the intermediate regions of 
many clusters, from where the entropy data are taken (see, 
e.g., Pratt et al. 2006, Pointecouteau et al. 2005). Thus, 
even though our comparison with the entropy data taken 
from Pratt et al. (2006) and from Ponman et al. (2003) 
can be quite safe, we more refined analyses of both the 
very inner parts (the cool cores at r ^ 0.1r2oo) and the 
external regions (those at r ^ r 2 oo) of clusters. The clus- 
ter cool cores that show quite high values of the magnetic 
field, sharp T jumps and probably not isentropic condi- 
tions (see, e.g., Piffaretti et al. 2005) and are also sites of 
complex AGN - IC gas interaction (see, e.g., Fabian 2005), 
are not considered in a complete self-consistent way here. 
Our approach to the description of the magnetized IC gas 
in this region needs a refined analysis of both the MVT 
and of the HE equation that is beyond the scope of this pa- 
per, and will be presented elsewhere. The external regions 
at r <L r2oo usually show a temperature decrease towards 
the outskirts of the cluster and thus, they also need a fur- 
ther analysis with respect to the isothermal model. The 
use of a polytropic model (i.e., an extension of our model 
that we have already worked out in the Appendix) will be 
also presented elsewhere. The study of the pressure and 
density structure of the cluster external regions is also rel- 
evant to set the appropriate boundary conditions for the 
IC gas and therefore it requires a more detailed analysis. 

Within the limits of its applicability, this first explo- 
ration indicates, nonetheless, that a more detailed phys- 
ical description of the IC gas properties that takes into 
account also the magnetic field might solve a number of 
problematic features in cluster structure and evolution. 
The predictions of our model can be subject to observa- 
tional tests. One of the most direct tests is to look for 
the amplitude of the magnetic field in clusters and groups 
through FR measurements that could directly prove or 
disprove the scaling of the magnetic field amplitude ~ 
j>o.5±o.i ^ na j- we p rec ii c t in order to recover the S — T and 
Lx — T relations. The few data pints available so far seem 
to be consistent with such a scaling (see also Cassano et 
al. 2006, for an independent argument). Detailed FR mea- 
sures in the inner ~ Mpc of groups and clusters will be 
available in the near future with the forthcoming SKA and 
LOFAR radio experiments. Some of the cluster will also 
likely show evidence for extended and diffuse radio emis- 
sion that will provide further information on the radial 
distribution of the magnetic field (see, e.g., Colafrancesco 
et al. 2005 for a more general discussion). Lastly, the com- 
bination of high-sensitivity X-ray and radio observations 
will be able to shed light on the physics of the intra-cluster 
gas and its associated magnetic field in both the cool cores 
and in the external regions of clusters at different redshifts, 
so that it will be possible to determine its detailed physical 
status as well as its cosmic evolution. 
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Appendix A: The hydrostatic equilibrium equation 
in the presence of B-field. 

The radial profile of the IC gas density p g can be obtained by 
solving the hydrostatic equilibrium equation once we specify 
the equation of state of the gas. In general one writes 



Vp = -p g V<j){r), 



(A.l) 



where p is the sum of all the pressure components in the cluster 
(IC gas, magnetic field, etc.) and <j>(r) is the gravitational po- 
tential. Under the assumption of spherical symmetry and local 
homogeneity of the gas, eq. (|A.l|) writes as 



1 dp 

p a dr 



d4 
dr 



GM(r) 



(A.2) 



where M(r) is the total mass enclosed in a sphere of radius r. 
The expression for M(r) depends on the assumed gravitational 
potential. In the /9-profile, one assumes that 

M(r) = 47r / 9 r s 3 {ln[x + (1 + x 2 ) 1 ' 2 ] - x(l + x 2 y 1/2 }. (A.3) 

Assuming a NFW universal profile, the mass M[r) writes as 



M{< r) = 47rp s (l + —*-) r l m (x) 

^ Mm 



(A.4) 



where m(x) describes the mass profile of the cluster. 

There are various choices for the gas equation of state: (a) 
the gas follows the Dark Matter, p g oc pdm, (b) the cluster 
is isothermal, T g =const, (c) the gas follows a polytropic law, 
T g oc pg -1 . In the following we will consider specifically the 
cases (b) and (c), and we will solve the hydrostatic equilibrium 
equation for the IC gas density with and without the presence 
of a magnetic field B(r). 

A.l. Isothermal case 

Let us consider first the case B = 0. 

In the case of a NFW mass density profile, eo, ()A.2[l can be 
written as: 



~Pgv) 



pm p 
Since 



GM(< r) 



dr 



M(< r) = M vir 



m(x) 



m(c) 

holds, it follows that 

dy g (x) _ GMvirpm p 
y g (x) k B T g r 2 

After integration we find 



m(x) 



i(c) 



x dy g (a 



o Va( x ) 

f x d y g {x) 



= -3 



= -3 



Tvir 



Tg m(c) 



dr. 



m(u) 



du, 



m(c) 



i(u) 



du. 



(A.5) 
(A.6) 

(A.7) 

(A.8) 
(A.9) 



from which 



P„(r) = p g (0)e ^ 



Jo 



obtains. The quantity T v i 



(A.10) 

is the virial tempera- 



ture and p 9 (0) is the central gas density. 

When B ^ 0, the hydrostatic equilibrium condition 
(eq |A.2[l is modified since an additional magnetic pressure term 
exists: 



d Pg (r,B) | dp B (r,B) 
dr 



GM(< r) 

~2 Poi.r,B), 



dr 
where 

B(r) ( xB t (pG)p a (r,Br 
and 

Pb oc B 2 . 

The HE condition in eq. ()A.ll|) can be written as: 



(A.ll) 
(A.12) 
(A.13) 



„ lr n s-x r dp 3 (r,B) 
p g {r,B) Gi + 

Or 87T 



Cl p 9 {r,B)-' 



d Ps (r,B) 
dr 



GM{< r) 



(A.14) 



and 



■p 9 (r,B) 



dp 3 (r,B) 2 
dr 



P^B)' 1 ^ 
GM(< r) 



2a 



2a - 1 



(A.15) 



where Ci = p g (0, B)/p g (0, B) oc k B T g and C 2 = B,(W 4 p g (z = 
0))~ a . It follows that 



yg r XiB) -x9y s ipBX + 



dr 



+ \2a-l) Cl Stt 



iCSS» e (r,B) 



dr 

After some algebra, the following equation 

2a \ C? 



dy g (x,B) + Pg (0,B) 

y s {x,B) 



p g (0,B) \2a-lJ 8tt 
p g (0,B) GM(< r) 



^ 2 -dy g (x,B) 2o! 1 



dr. 



Pg (Q,B) r 2 
is found, and after integration one obtains: 



(A.17) 



lny fl (r,B) + rj'{y g {r, B) 2 



D = 



\ny g (r,0) 



M^(B)2 47r r^ r 



G M 3 " 



" Pext 



C 2 ( 2a \ R x2a-l prn 

= ^{2*—l) P3i0 > B) l^Tg 



(A.18) 



(A.19) 



where y g (r, 0) is defined by equation eo. (|A,10[l . The solution 
of this equation provides the radial profile of the gas density 
once the constants p a (0, B) and T g (B) are known. The central 
gas density is set as 



p g (0,B)=p g (0,0), 
while the temperature 

T g (B)=T g (0) + 



M 2 . 

v%r 



G Ml, 



(A.20) 



(A.21) 
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is obtained from the MVT (see Colafrancesco & Giordano 
2006a for a derivation). For a given cluster gas mass, the den- 
sity profile should be re-normalized as 



Mg }V i r (B) = Mg yV i r (B = 0) 

from which the condition 

/ dxx 2 y g (x,0) 



p 9 (B,0)=p 9 (0,0) 



obtains. 



/ dxx 2 y g (x,B) 



(A.22) 



(A.23) 



A. 2. Poly tropic case 

Let us consider first the case B = 0. 
Assuming 



Pa <x Pg, 



eq. (|A.2[) writes as: 



KiPair)- 1 



dpgjr) 
dr 



GM(< r) 



(A.25) 



with K\ = p 9 (0)/p 9 (0) 7 . From this one obtains 
Ki- 



7-1 
and 

dr 



7 rc^ dfrfr-r- 1 _ GM(< r) 



dr 



7- 1 \ GA/(< r)pm p 
7 J fc B T g (0)r2 ■ 



(A.27) 



From eg . [|A.6|I and after integration between the limits r — 
e r, one finds the profile 



Pg( r ) = p g (0) 



1-3 



T vir ( 7 — 1 



r g (o) V 7 / m(c) y w : 



-du 



with 

Tvir 



3 c i&r 1 



7 m(c) 
m(c) 



+ (7-l)l&| 



m(u) 



(/it 



7 = 1.15 +0.01(c- 6.5), 



(A.28) 



(A.29) 



(A.30) 



and S* = dlnydm(x) /dlnx\ c . 

For B^O, ea. (IA.ll[) is substituted by the following equa- 
tion: 

K iP a (r,B) ggU + ^p 9 (r,B) 



a^vB) 

dr 



and 



7 \ dp g (r,B) 



7-1 



GM(< r) 



(A.31) 



7 — 1 J dr 
dp g (r,B) 2 -< 



dr 



= -Ki 



+ Al Stt Ua-1 
4 GM(< r) 



(A.32) 



with _R"2 = C2. Then, it follows 

3y s (r,5) 7 ^ , /^7-X\/ 2a \^|p 9 (0,B) 2 



(9r 



2a -1/ 8tt p 9 (0,B) 



9y 9 (r,B) 



7-l\ P 9 (0,B) GM(< r) 



dr \ 7 J p 9 (0, B) r 2 

After integration one finds 



(A.33) 



y]- 1 (x,B)+r t y 2 g a - l (x,B) = (1 + V )y]- 1 (x,B = 0), (A.34) 
where 



Pi3(r = 0,B) 2a 7-I 
p 9 (r = 0,5) 2a -1 7 ' 



(A.35) 



To also solve numerically eq. ()A.34|l we have to know the quan- 
tities 7, T 9 (0, B) and p g (0,B). The quantity 7 is independent 
of the magnetic field B\ the quantities T g (0,B) and p g (0,B) 
are related to the corresponding quantities for B — by the 
relations: 



(A.24) T g (0,B) 



T 9 (0,0) 



I + 77 ' 

As a consequence, it follows that: 
p 9 (0,0) 



pg(0,B) 



(l + V) 1 



(A.36) 



(A.37) 
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